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Dynamical bar instability in a relativistic rotational core collapse 



in ; 
o . 
o ■ 

(N ■ 

a ■ 



Motoyuki SaijcQ 

Department of Physics, Kyoto University, Kyoto 606-8502, Japan 
(Received 28 February 2005; Published 31 May 2005) 

We investigate the rotational core collapse of a rapidly rotating relativistic star by means of a 3+1 
hydrodynamical simulations in conformally flat spacetime of general relativity. We concentrate our 
investigation to the bounce of the rotational core collapse, since potentially most of the gravitational 
waves from it are radiated around the core bounce. The dynamics of the star is started from a 
differentially rotating equilibrium star of T/W ~ 0.16 (T is the rotational kinetic energy and W 
is the gravitational binding energy of the equilibrium star), depleting the pressure to initiate the 
collapse and to exceed the threshold of dynamical bar instability. Our finding is that the collapsing 
star potentially forms a bar when the star has a toroidal structure due to the redistribution of the 
angular momentum at the core bounce. On the other hand, the collapsing star weakly forms a bar 
when the star has a spheroidal structure. We also find that the bar structure of the star is destroyed 
when the torus is destroyed in the rotational core collapse. Since the collapse of a toroidal star 
potentially forms a bar, it can be a promising source of gravitational waves which will be detected 
in advanced LIGO. 

PACS numbers: 04.25.Dm,04.30.Db,04.40.Dg,97.60.-s 
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I. INTRODUCTION 

Rotational core collapse occurs in the first phase of the 
collapse-driven supernova explosion. The collapse takes 
place in a dynamical time scale and the bounce occurs 
when the central part of the star collapses into a nuclear 
density and then stiffens its equation of state. The final 
outcome of the collapse forms a compact object, such as 
a neutron star or a black hole, which is a candidate for 
sources of gravitational waves 0. There are two differ- 
ent approaches to investigate the rotational core collapse. 
One is to take into account the realistic picture of the ro- 
tational collapse such as realistic equation of state of the 
neutron star and/or neutrino effect in Newtonian gravity 
[e.g.0,UQ. The other is to illustrate simplified physics 
in relativistic gravitation [e.g.0,@- Since one of our aims 
in this paper is to investigate the possibility of gravita- 
tional wave sources in rotational core collapse, we should 
at least take into account relativistic gravitation, which 
significantly affects the quantitative behavior of gravita- 
tional waves from rotational core collapse 0. 

Dynamical bar instability in a rotating equilibrium 
star takes place when the ratio f3(= T/W) between ro- 
tational kinetic energy T and the gravitational bind- 
ing energy W exceeds the critical value /3d yn - Deter- 
mining the onset of the dynamical bar-mode instabil- 
ity, as well as the subsequent evolution of an unstable 
star, requires a fully nonlinear hydrodynamic simula- 
tion. Simulations performed in Newtonian gravity [e.g. 
HUEmEBElHEIIl have shown that (3 dyn 
depends only very weakly on the stiffness of the equa- 
tion of state. /?dyn becomes small for stars with high 
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degree of differential rotation 0, 0, |2(| • Simulations in 
relativistic gravitation 0, [2^ have shown that Pdyn de- 
creases with the compaction of the star, indicating that 
relativistic gravitation enhances the bar-mode instability. 
The dynamical bar instability potentially occurs during 
the collapse since the nondimensional value /? scales as 
in the dimensional analysis where R is the radius of 
the star. Let us briefly discuss a picture of core bounce. 
When the star collapses, (3 increases and the star starts 
to deform its shape to form a bar when the star exceeds 
the critical value of dynamical instability. During the 
bounce phase, (3 falls below the threshold of dynamical 
instability, and the star cannot deform furthermore. In 
such a situation, what happens to the nonaxisymmetric 
deformation of the star? 

There are several papers that investigate rotational 
core collapse as a potential source of gravitational waves 
in axisymmetric spacetime in Newtonian gravity 0, l23| , 
in conformally flat spacetime |5|,|6|,|24j, an d in full general 
relativity [23, [2a> l27l to estimate the amount of gravi- 
tational radiation [28j |. Recently, 3D calculations in full 
general relativity have been established to investigate the 
nonaxisymmetric deformation of the star in rotational 
core collapse. Duez, Shapiro, and Yo (2^| investigated 
the collapse of a differentially rotating n = 1 polytropic 
star in 3D by depleting the pressure and found that the 
collapsing star forms a torus which fragments into non- 
axisymmetric clumps. Shibata and Sekiguchi [3H | inves- 
tigated in rotational core collapse and found that a burst 
type of gravitational waves was emitted. In addition, 
they argued that a very limited window for the rotating 
star satisfies to exceed the threshold of dynamical insta- 
bility in the core collapse. Zink et al. [3l| presented a 
fragmentation of an n = 3 toroidal polytropic star to a 
binary system inducing m = 2 density perturbation and 
claimed that they found a new scenario to form a binary 
black hole. 

Our purpose in this paper is the following twofold. One 
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is to investigate a necessary condition and mechanism to 
enhance the dynamical bar instability in a collapsing star. 
Brown performed long time integration of the rota- 
tional core collapse to investigate whether a dynamical 
bar can significantly form during the evolution. He found 
that the dynamical bar instability sets in at /3d yn ~ 0.23, 
which is far below the standard value /?dyn ~ 0.27. He 
found that the role of dynamical bar in rotational core 
collapse may be quite different from that in the equilib- 
rium star. He also found that the bar instability grows 
slowly in core bounce. His finding can be interpreted that 
bar instability is initiated by the interaction between the 
core part and the surrounding part of the star. Therefore 
the "dynamical" bar instability in the collapsing star is 
quite different from that in the equilibrium state. 

The other is the importance of probing whether the ro- 
tational core collapse becomes a promising source of grav- 
itational waves. Direct detection of gravitational waves 
by ground based and space based interferometers is of 
great importance in general relativity in astrophysics, 
and in cosmology. Once a bar has formed in the neutron 
star, we may expect quasi-periodic gravitational waves 
in the kHz band, which may be detectable in advanced 
LIGO [[[. Rampp, Miiller, and Ruffert (3j| investigated 
the rotational core collapse using parametric equation of 
state in Newtonian gravity, including the thermal pres- 
sure and the nuclear matter when the density exceeds 
the threshold of nuclear density. They found that al- 
though the rotating star has a nonaxisymmetric defor- 
mation, the amplitude of gravitational waves does not 
significantly grow during the deformation. In fact, they 
compared the gravitational waveform computed in their 
simulation with that in the previous 2D calculation j2^| 
and found that it is quite similar to each other. There- 
fore, even if the nonaxisymmetric instability takes place 
during the collapse, the bar does not significantly grow 
at the core bounce. What else do we need to enhance the 
bar formation in the rotational core collapse? 

Our interest is to focus on the core bounce of the ro- 
tational core collapse in order to investigate the angular 
momentum redistribution and to investigate the possi- 
bility of gravitational wave sources. Accordingly we con- 
struct a mimic model to investigate the above two issues. 
We deplete the pressure of the equilibrium star to initiate 
collapse and bounce. We also concentrate on the struc- 
ture of the star, spheroidal and toroidal, to investigate 
the dynamical bar instability in the collapsing star. Stel- 
lar collapses and mergers may also lead to differentially 
rotating stars [e.g. For the coalescence of binary 

irrotational neutron stars [3^ |3(| |23, the presence of 
differential rotation may temporarily stabilize the "hy- 
permassive" remnant, which constructs a toroidal struc- 
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ture and may therefore have important dynamical effects. 
Although they use stiff equation of state (n — 1) it is pos- 
sible to construct a toroidal star in nature. 

This paper is organized as follows. In Sec. ^we sum- 
marize our basic equation of relativistic hydrodynamics 
in conformally flat spacetime. In Sec. IIIII we show our 
numerical results of dynamical bar instability in rota- 
tional core collapse, and summarize our findings in Sec. 
IIVI Throughout this paper we use geometrized units 
(G = c = 1) and adopt Cartesian coordinates (x, y, z) 
with the coordinate time t. Greek and Latin indices run 
over {t,x,y,z) and (x,y,z), respectively. 



II. RELATIVISTIC HYDRODYNAMICS IN 
CONFORMALLY FLAT SPACETIME 

In this section, we describe the basic equations in con- 
formally flat spacetime [e.g. |H, |3^, 01 ■ We solve the 
fully relativistic equations of hydrodynamics, but neglect 
nondiagonal spatial metric components. 



A. The gravitational field equations 

We define the spatial projection tensor ^ v = g^ v + 
n^n 1 ' , where is the spacetime metric, = 

(I/a, — (3 l / a) the unit normal to a spatial hypersurface, 
and where a and (3 l are the lapse and shift. Within a 
first post-Newtonian approximation, the spatial metric 
9ij = lij may always be chosen to be conformally flat 

7y=V> 4 % (2.1) 

where ip is the conformal factor [see 0, ^3 • The space- 
time line element then reduces to 

ds 2 = {-a 2 + PkP k )dt 2 + 2/3 i dx i dt + ip 4 6 ij dx i dx j . 

(2.2) 

We adopt maximal slicing, for which the trace of the 
extrinsic curvature vanishes, 

K = i l ^K l3 = 0. (2.3) 

The gravitational field equations in conformally flat 
spacetime for the five unknown a, (3 l , and tp can then 
be derived conveniently from the 3+1 formalism. The 
equation for the lapse a, shift /3 1 , and conformal factor ip 
with a maximal slicing implies dtK — 0, shall be written 
as 



A(«V) = 2iraiP 5 { Pli + 2S) + -aip 5 KijK i3 , (2.4) 
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S u A(3 l + ^dtdifi 1 = WiraJt + (dj In 
A-0 = -2ttiP 5 Ph - -r r 7\, ; /\'<. 



where 5 — 7jfeT J ' , A = S^didj is the flat space Lapla- 
cian and J t = —r\} lr f^T lhV is the momentum density. 
In the definition of Ji, T^ v is the stress energy tensor, 
Ph = n^n^T^ is the mass-energy density measured by a 
normal observer. 



B. The matter equations 

For a perfect fluid, the energy momentum tensor takes 
the form 

T" u = p (l + e + ^ u"u" + Pg" u , (2.7) 

where p is the rest-mass density, e the specific internal 
energy, P the pressure, and w M the four-velocity. We 



(^) ) ( dlPJ + 5il5jkdk P l ~ ^ W) . ( 2 - 5 ) 

(2.6) 

I 

adopt a T-law equation of state in the form 

P=(T-l)pe, (2.8) 

where T is the adiabatic index which we set T = 5/3, 3/2, 
7/5 in this paper. In the absence of thermal dissipation, 
Eq. H2.8JI . together with the first law of thermodynamics, 
implies a polytropic equation of state 

P = K p 1+1/n , (2.9) 

where n — 1/(T — 1) is the polytropic index and k is a 
constant. 

From V ll T pv — together with the equation of state 
(Eq. 12.81 ). we can derive the energy and Euler equations 
according to 



| + M.^4 W v )l ( , 10) 

~, + a ■ = -oeip b (P + P v i S ),i - ptoafaj + p*u 3 p J i + r (2.H) 



where 



pav^i/j 6 , 

(i + r e y, 

(1 +re)w 2 , 



(2.12) 

(2.13) 

(2.14) 
(2.15) 
(2.16) 



and u M , P V i S is the 3- velocity, 4- velocity, pressure vis- 
cosity, respectively. Note that we treat the matter fully 
relativistically; the conformally flat approximation only 
enters through simplifications in the coupling to the grav- 
itational fields. As a consequence to treat shocks we also 
need to solve the continuity equation 



<9/9» d(p*v l ) 



dt 



dx l 



0. 



(2.17) 



separately. 



C. Numerical techniques for solving gravitational 
field equations 

We have reduced Einstein equations in conformally flat 
spacetime to ten elliptic equations for ten variables (Bj, 
X, ip, aM),Pi, rf) using the same techniques in the previous 
paper [4(j 



Ax 



STT^JiX 1 EE 47r5 x , 

-2tt^ 5 ph - \^- 7 A t3 A^ ee 4^5^,(2.20) 



(2.18) 
(2.19) 



AV> 

A(o^) = 2natp(p K + 25) + ^aip~ 7 A tJ A i] 

EE 47r5, 
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APi = iiraJi ee 47r5p s , 
Ar] = —AiraJiX 1 = 4^5^. 



(2.21) 
(2.22) 
(2.23) 



We use the asymptotic fall off behavior for metric 
quantities at large radius in order to set an appropriate 
boundary condition at the grid edge [40| . The definition 
of P>i, x, Pi, V and the basic procedure to solve ten elliptic 
equations are given in Ref. [40j . Since we parallelized our 
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TABLE I: Relativistic Rotating Equilibrium stars 



Model 


n a 


A b 


Rp 1 Re C 


niax d 

Po 




Rc e 


M f 


J 9 




T/W h 




Rc/M 


p i 

* dcp 


I-a 


1.5 


0.5 


0.467 


8.26 x 10" 


-4 


6.56 


1.64 x 10" 1 


4.49 x 10" 


2 


1.59 x 10" 


i 


3.99 x 10 1 


70% 


I-b 


1.5 


0.8 


0.467 


1.36 x 10" 


-3 


7.22 


1.79 x 10" 1 


5.25 x 10" 


-2 


1.59 x 10" 


i 


4.04 x 10 1 


70% 


Il-a 


2.0 


0.5 


0.450 


1.10 x 10" 


-4 


2.53 x 10 1 


6.29 x 10 _1 


6.33 x 10" 


1 


1.63 x 10" 


i 


4.03 x 10 1 


60% 


Il-b 


2.0 


0.6 


0.417 


2.19 x 10" 


-4 


2.64 x 10 1 


6.62 x 10" 1 


6.69 x 10" 


1 


1.62 x 10" 


i 


3.98 x 10 1 


60% 


Ill-a 


2.5 


0.3 


0.417 


1.05 x 10" 


-S 


8.63 x 10 1 


2.16 


6.90 




1.59 x 10" 


i 


3.99 x 10 1 


50% 


Ill-b 


2.5 


0.42 


0.417 


2.24 x 10" 


-5 


8.54 x 10 1 


2.13 


6.58 




1.60 x 10" 


i 


4.00 x 10 1 


50% 



"Polytropic index 

'Parameter of the degree of differential rotation 

c Ratio of the polar proper radius to the equatorial proper radius 

''Maximum rest-mass density 

e Equatorial circumferential radius 

^Gravitational mass 

9 Angular momentum 

h T: Rotational kinetic energy; W: Gravitational binding energy 
'Ratio of pressure depletion to the equilibrium star 



code, we change the method to solve the elliptic equation 
to PCG method [e.g. E|. 

We monitor the rest mass Mo, gravitational mass M, 
and the angular momentum J 



M 
M 



= I p*d 3 x, 



(2.24) 



[(p + pe + P)(au t ) 2 - P] i/; 5 



1 

16tt 



d 3 x, 



J = 



(2.25) 



(2.26) 



during the evolution. We also compute rotational kinetic 
energy T, proper mass M p , gravitational binding energy 
W as 



(xJ y - yJ x )'ilj 6 d 3 x, 



Mr, 



T 



pu t (l + e)^gd 3 x 
p*(l+e)d 3 x, 



- I nriv^gd'x 



fl(xj y — y.J x )ip d x 



W = M p + T-M, 



(2.27) 

(2.28) 
(2.29) 



where 17 the angular velocity of the star. 

Since we use a polytropic equation of state at t = 0, 
it is convenient to rescale all quantities with respect to 



k. Since k™/ 2 has dimensions of length, wc introduce the 
following nondimensional variables [e.g. |44| 



-n/2„ 



-™/2, 



t = K n < t, x — k "i"x, y = K "'"y, 
z = K- n ' 2 z, Q = K n / 2 fl, M = k-"/ 2 M, (2.30) 
R = K ~ n/2 R, J = n- n J. 



Henceforth, we adopt nondimensional quantities, but 
omit the bars for convenience (equivalently, we set n — 

Our code is based on the conformally flat hydrodynam- 
ics scheme of Ref. 40], to which the reader is referred 
for a more detailed description, discussion and test. Wc 
choose the axis of rotation to align with the z axis, and 
assume planar symmetry across the equator. The equa- 
tions of hydrodynamics are then solved on a uniform grid 
of size 200 x 200 x 60. We terminate our simulations after 
a sufficient number of central rotation periods (between 8 
and 16) in order for us to detect dynamical instabilities. 
Because of our flux-conserving difference scheme the rest 
mass Mo is conserved up to a round-off error, except if 
matter leaves the computational grid (which was never 
more than 0.01% of the total rest mass). In all cases re- 
ported in Sec. lIIII the total gravitational mass M and the 
angular momentum J were conserved up to ~ 0.1% and 
less than about 5% of their initial values. 



III. ROTATIONAL CORE COLLAPSE 

We basically follow the scheme of Ref. to con- 

struct a differentially rotating equilibrium star. The de- 
tailedprocedure we used for constructing it is given in 
Ref. I43. We construct the equilibrium star under the 
condition of fixing (3 and R c /M, where R c is the circum- 
ferential radius. This is because we keep approximately 
the same condition to extend the threshold of dynamical 
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FIG. 1: Evolution of quadmpole diagnostics in 6 different collapsing stars. Roman character in the panel corresponds to the 
model type in Table [I] Solid lines and dash lines represent toroidal and spheroidal stars, which correspond to model a and b in 
Table respectively. 
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FIG. 2: Snapshot of the coordinate density along the x-axis for 6 different stars (see Table QJ. Solid lines and dash lines 
represent t/P c = 1(a) (9.61, 16.0), 1(b) (9.67, 15.3), 11(a) (2.39, 12.0), 11(b) (3.25, 12.0), 111(a) (3.19, 7.77), 111(b) (3.20, 7.77), 
respectively. 
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FIG. 3: Density contour in the equatorial plane of 2 collapsing stars (Model I [a] and [b] of Table 0. Snapshots are plotted 
at the parameter (t/t dyn , p* msx ) = I(a)-i (9.61, 9.13 x 10" 3 ), I(a)-ii (16.0, 6.04 x 10" 3 ), I(b)-i (9.67, 1.33 x 10" 2 ), I(b)-ii (15.3, 
1.57 x 10~ 2 ), respectively. The contour lines denote coordinate densities p* — p„ ax x 10~ 0,220 ' 16-1 ' (i = !,-■■, 15). 



bar instability by depleting the pressure of the star to ini- 
tiate collapse. We set one toroidal star (corresponds to 
model a in Table [IJ and one spheroidal star (corresponds 
to model b in Table P) for each polytropic index to focus 
on the structure of the star. 

We choose the rotation profile of the star as |45j 



u u u 



A 2 (n c -n), 



which represents in the Newtonian limit (it* 
ru 2 il), so-called j-constant law, as 



Q = 



A 2 n c 



A*' 



(3.1) 

1, u v — > 

(3.2) 



where A is a parameter representing the degree of dif- 
ferential rotation, w is the cylindrical radius of the star. 
Since A has a dimension of length, we normalize it with 
the proper equatorial radius R e , (A = R e A). We summa- 
rize our 6 different initial data sets of relativistic rotating 
equilibrium stars in Table [i] 

To monitor the development of m — 2 modes we com- 



pute a "quadrupole diagnostic" j4y| 

1 f (x 2 - y 2 ) + i(2xy) 3 



r 



X, 



where a bracket denotes the density weighted average. In 
the following we only plot the real parts of Q. 

To enhance any m = 2 instability, we disturb the initial 
equilibrium density p cq by a nonaxisymmetric perturba- 
tion: 



p = Peq I 1 + 5 



y 



R 2 



(3.4) 



with 8 = 0.01 in all our simulations. 

As for computing the gravitational waveform, we use 
the same method that we used in the previous paper |22j . 
For observers along the rotational axis (z-axis), we have 



rft_i_ 1 d 2 

— - — — (/ -J ) 

M ~ 2Mdt 2{ xx yvh 



rhy, 
M 



Mdt 2 xv ' 



(3.5) 
(3.6) 
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FIG. 4: Density contour in the equatorial plane of 2 collapsing stars (Model II [a] and [b] of Table Pi. Snapshots are plotted at 
the parameter (t/t dyn , p* msx ) = II(a)-i (3.19, 1.13 x 10~ 3 ), II(a)-ii (12.0, 4.62 x 1(T 3 ), II(b)-i (3.25, 3.39 x 1(T 3 ), II(b)-ii (12.0, 
8.71 x 10~ 3 ), respectively. The contour lines denote coordinate densities p" — p^, ax x io~ 0,200 ' 16 ~'' (i = !,■■■, 15). 



where h + and h x are the two polarization modes of grav- 
itational waves, r is the distance to the source, Zy the ap- 
proximate quadrupole moment of the mass distribution 
defined as 

Iij = J p*x z x j d 3 x. (3.7) 

We show the quadrupole diagnostic Q throughout the 
evolution in Fig. ^ We determine that the system is 
stable to m = 2 mode when the quadrupole diagnostic 
remains small throughout the evolution. We also deter- 
mine that the system is unstable when the diagnostic 
grows exponentially at its early stage of evolution. It is 
clearly shown in Fig. 2] that the star is more unstable to 
the bar mode when the structure of the star is toroidal at 
its equilibrium stage than when it is spheroidal. Even in 
a spheroidal case, the dynamical instability sets in after 
the core bounce in case of model 1(b) (see Table 01. Note 
also that the diagnostic damps out especially for the soft 
equation of state (n = 2, 2.5). The above effect corre- 
sponds to the destruction of the toroidal structure of the 
star (see Fig. 



We show the density contour in the intermediate stage 
and in the final stage of the rotational core collapse in 
Figs. |21 El In the intermediate stage of the rotational 
core collapse, the bar instability grows due to the nonax- 
isymmetric instability. The deformation rate of the bar 
is significantly higher for the case of a toroidal star at 
equilibrium than that of a spheroidal star. At the termi- 
nation of integration, the equatorial shape of the star is 
almost spherical except for the case of model 1(a). The 
coincidence of the destruction of the toroidal structure 
in Fig. [5] at the termination of integration indicates that 
the structure of torus plays a significant role in enhancing 
bar instability. 

We show the gravitational waveforms along the ob- 
server in the rotational axis in Fig. The gravitational 
waves amplify during the bar formation. Since the dy- 
namical bar plays an important role in our case of ro- 
tational core collapse, our result is quite different from 
the previous picture of the 3D Newtonian results [3j| • I n 
fact, the amplitude in our calculation observed from the 
rotational axis grows about 20 times larger than that of 
core bounce which has a peak in the waveform. Note 
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FIG. 5: Density contour in the equatorial plane of 2 collapsing stars (Model III[a] and [b] of Table||J. Snapshots are plotted 
at the parameter (t/t dyn , p^) = III(a)-i (3.19, 2.70 x 1(T 4 ), III(a)-ii (7.77, 4.02 x 10~ 3 ), III(b)-i (3.20, 3.28 x 10" 4 ), III(b)-ii 
(7.77, 1.92 x 10~ 3 ), respectively. The contour lines denote coordinate densities p" = p^ax x 10~°- 220(16_l) (i = !,-■■ , 15). 



also that the quasiperiodic waves retain several rotation 
periods for stiff equation of state, due to the persistence 
of the toroidal structure. 

We also show (3 in the evolution in Fig. which is 
regarded as a diagnostic of the dynamical bar instability 
in the equilibrium star. Although (3 behaves quite sim- 
ilar in the different polytropic index, the bar structure 
persists for at least several rotation periods in case of the 
model 1(a), which corresponds to the persistence of the 
toroidal structure (see Fig 121) • We also show the distri- 
bution of the angular velocity in the intermediate stage 
and at the termination of the integration in Fig. |H1 A 
sharp dip at the center in the angular velocity of model 
1(a) at t = 9.61P C potentially has a redistribution of the 
angular momentum. 



IV. DISCUSSION 

We investigate the role of dynamical bar instability 
in rotational core collapse by means of hydrodynamic 
simulations in conformally flat approximation in general 



relativity. We specifically focus on the structure of the 
star to see whether it takes a significant role in enhancing 
dynamical bar instability. 

We find that the structure of the star takes a signif- 
icant role in enhancing the dynamical bar instability at 
core bounce. Since the angular velocity of the collapsing 
star has already reached the maximum inside a certain 
cylindrical radius to produce a toroidal structure, the an- 
gular momentum can only shift outward at the bounce. 
For a spheroidal star, the angular momentum can shift 
both inward and outward at bounce since it does not 
reach the "Keplarian" . This means that rotational core 
collapse for the spheroidal case cannot significantly break 
the central core of the star. Consequently, in case of a 
toroidal star a bar structure is easily constructed dur- 
ing the evolution. Note that for a soft equation of state 
(n = 2.0, 2.5) the amplitude of gravitational waves de- 
creases, as the torus is destroyed in the rotational core 
collapse. Therefore the torus is the key issue to trigger 
the bar formation. 

Once a bar has formed, the amplitude of gravitational 
waves has significantly increased due to the nonaxisym- 
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FIG. 6: Gravitational waveforms from a distant observer along the rotational axis of the equilibrium star. Roman character 
in the panel corresponds to the model type in TableQ Solid lines and dash lines represent toroidal and spheroidal stars, which 
correspond to a and b in Table HJ respectively. 
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metric deformation of the star. The previous 2D calcula- 
tion shows that a peak amplitude of gravitational waves 
comes from the core bounce of the star, and its behav- 
ior coincides with that of 3D calculation. In our results, 
gravitational radiation is dominantly generated in the bar 
formation process. Our results of the soft equation of 
state qualitatively agree with the results in full general 
relativity 30]. Thus, the characteristic amplitude and 
frequency of gravitational waves in the collapsing star 
can be written as 



/gw 



-J- = 100 [Hz] ( A M 



3/2 



, (4-1) 



4.8 x 10" 



-23 ( M \ /^lOMpc^ 



rh/M\ 
0.01 ) 



(4.2) 



Therefore, gravitational waves from a dynamical bar in 
a collapsing neutron star will be detected in advanced 
LIGO. 

Finally we mention the validity of the conformally flat 
approximation in our 6 different collapsing stars. The 
approximation has two issues: that it contains the first 
post-Newtonian order of general relativity, and that it 
neglects gravitational radiation. First, we investigate the 
central lapse a c at core bounce to check whether the first 
post-Newtonian order of general relativity is reasonable 
approximation in our model. The central lapse can be 
expanded by the compaction (characteristic mass M„ r 
and radius R SIC ) of the source [Eq. (19.13) of Ref. |47| ] 



TABLE II: Magnitude of higher order correction for the cen- 
tral lapse at bounce 



Model t bco /P c J 


a c k 




(M src /i? src ) 3 [%] 


I-a 


1.22 


0.872 


0.137 


0.257 


I-b 


1.00 


0.848 


0.166 


0.457 


Il-a 


1.05 


0.843 


0.172 


0.509 


Il-b 


0.88 


0.806 


0.218 


1.04 


Ill-a 


1.08 


0.835 


0.181 


0.593 


IU-b 


0.90 


0.803 


0.222 


1.09 



3 tbcc : Bounce time 
^Central lapse 

; M arc : Characteristic mass of the source; -R arc : 
radius of the source 



Characteristic 



as 



otc. = 1 - 



M S1 



1 (M* 

2 V R SI 



\Rs 



(4.3) 



Note that the shift appears in the gtt component of the 
4-metric at the second post-Newtonian order of general 
relativity. All models of our 6 collapsing stars bounce at 
a c > 0.8, and therefore the value of central lapse corre- 
sponds to the compaction of the source M src /i? src < 0.22 
from Eq. (|4.3|l . For each model, we summarize the 
magnitude of higher order correction [(i? S rc/-^src) 3 term] 
from the first post-Newtonian order in Table [H] Note 
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that the central lapse at bounce is the minimal one 
throughout the evolution. Since the magnitude of the 
higher order correction is (M src /i? src ) 3 < 0.011, we can 
roughly claim that the first post-Newtonian gravity is 
a satisfactory approximation to describe the full general 
relativity in the relative error rate of few percent in our 
calculation, depending on the coefficient of the higher 
order term. Second, it is also a quite satisfactory ap- 
proximation to neglect gravitational radiation in our dy- 
namics, since our main target in this paper is to focus on 
the enhancement of the dynamical bar instabilities which 
occur in dynamical time scale. Gravitational waves affect 
the whole dynamics in secular time scale, which is much 
longer than the dynamical time scale, and hence we can 
safely neglect such effect in this paper. Therefore, con- 
formally flat approximation is a satisfactory one in our 
computation. 
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